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^ , Abstract - Analytical formulae for de Haas- van Alphen (dHvA) oscillations in linear chain of cou- 

O ■ pled two-dimensional (2D) orbits (Pippard's model) are derived systematically taking into account 

the chemical potential oscillations in magnetic field. Although corrective terms are observed, basic 
[a) and magnetic breakdown-induced {(5 and 2(5 - a) orbits can be accounted for by the Lifshits- 
^ ■ Kosevich (LK) and Falicov-Stachowiak semiclassical models in the explored field and temperature 

\^ I ranges. In contrast, the 'forbidden orbit' p - a amplitude is described by a non-LK equation in- 

^\ , volving a product of two classical orbit amplitudes. Furthermore, strongly non-monotonic field and 

\Q ■ temperature dependence may be observed for the second harmonics of basic frequencies such as 

CO \ 2a and the magnetic breakdown orbit /3 + a, depending on the value of the spin damping factors. 

These features are in agreement with the dHvA oscillation spectra of the strongly 2D organic metal 
6'-(ET)4CoBr4(C6H4Cl2). 



Introduction. — Quasi- two-dimensional (q-2D) multiband organic compounds (ET)2X| 
(where ET stands for the bis-ethylenedithio-tetrathiafulvalene molecule and X is a monova- 
lent anion) exhibit very rich phase diagrams with ground states ranging from antiferromag- 
netic insulator to Fermi liquid. In particular, superconducting and magnetic ordered phases 
are observed close to each other which is reminiscent of cuprates and iron-based supercon- 
ductors (for a recent review, see [1 and references therein). These compounds share the 
same Fermi surface (FS) topology which is an illustration of the linear chain of coupled orbits 
introduced by Pippard in the early sixties [2 . Such FS is liable to give rise to a network of 
orbits coupled by magnetic breakdown (MB) in large magnetic fields. Quantum oscillations 
in such systems have been extensively studied during the two last decades (for a review, see 
e.g. [3] and references therein). The main feature of the oscillation spectra of these com- 
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pounds is the presence of frequency combinations that are forbidden in the framework of the 
semiclassical model of Fahcov-Stachowiak [4l. Theoretical studies have demonstrated that 
these Fourier components result from both the formation of Landau bands due to coherent 
MB, instead of discrete levels, and the oscillation of the chemical potential enabled by the 
2D character of the FS. However, easy to handle analytic tools necessary to quantitatively 
account for the data are still lacking. As a matter of fact, comprehensive theoretical cal- 
culations taking into account coherent MB and oscillation of the chemical potential have 
been reported in Ref. [5] and implemented in de Haas- van Alphen (dHvA) data of the above 
mentioned compound. Nevertheless, as pointed out in Ref. [5 , equations are complex and 
data could be analyzed only numerically. 

The aim of this paper is therefore to provide analytical formulae accounting for quantum 
oscillation spectra of these systems. In the first step, analytical calculation of the Fourier 
amplitude of the various dHvA oscillations predicted for the FS illustrating the Pippard's 
model is reported. In the second step, the FS of the recently synthesized strongly 2D charge 
transfer salt ^-(ET)4CoBr4(C6H4Cl2) [6] is presented and, in the third step, analysis of the 
observed field- and temperature-dependent quantum oscillation spectra are reported. An 
excellent agreement with the calculations is obtained. 

Experimental. — The studied crystals were synthesized by electrocrystallization tech- 
nique as reported in Ref. [GJ. They are elongated hexagonal platelets with approximate 
dimensions 1.5 x 0.5 x 0.04 mm^. X-ray diffraction data were collected both at room tem- 
perature in Chernogolovka with a KM-4 single-crystal diffractometer (Kuma Diffraction) 
and at 100 K with an Xcalibur diffractometer (Oxford Diffraction) at the Labor atoire de 
Chimie de Coordination in Toulouse. Tight-binding band structure calculations were based 
upon the effective one-electron Hamiltonian of the extended Hiickel method [7] , as reported 
e.g. in Ref. [8|. Magnetic torque and magnetoresistance were measured in pulsed magnetic 
fields of up to 55 T with a pulse decay duration of 0.32 s. Magnetic torque measurements 
were performed with a commercial piezoresistive microcantilever [9] in the temperature 
range from 1.4 K to 4.2 K. The size of the sample studied was approximately 0.12 x 0.1 
X 0.04 mm^. Variations of the cantilever piezoresistance was measured with a Wheatstone 
bridge with an ac excitation at a frequency of 63 kHz. The angle between the normal to 
the conducting plane and the magnetic field direction was = 7°. Magnetoresistance was 
measured thanks to a contactless radio frequency measurement technique based on tunnel 
diode oscillator as reported in Ref. A pair of counter- wound coils made with copper 

wire of 50 jiia in diameter was used. As pointed out in [11], in the case where the inter- 
layer resistivity is much larger than the in-plane resistivity, the contribution of the former 
parameter to the measured signal is negligible. The time evolution of the oscillator fre- 
quency during the field pulse is deduced from successive short-term fast Fourier transforms 
of the raw signal. In contrast, analysis of the oscillatory magnetoresistance and magnetic 
torque is based on discrete Fourier transforms of the field-dependent data, calculated with 
a Blackman window after removal of a non-oscillating background. Zero-field interlayer and 
in-plane 4-point resistance was measured with 1 jiA ac current at 77 Hz. Electrical contacts 
were made using annealed platinum wires of 20 /im in diameter glued with graphite paste. 

Calculation of de Haas- van Alphen amplitudes. — We consider a two-band sys- 
tem with band extrema Ao(i) (see Fig. [He)) and effective masses mo(i). In the low temper- 
ature limit, the grand potential O of a two-dimensional slab with area A is given by a series 
of harmonics p of classical frequencies : 



__f^(T,;,) = -_(^-Ao) --(M-Ai) 

p>l r; 
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Fig. 1: (Colour on-line) (a) Crystal structure of ^-(ET)4CoBr4(C6H4Cl2) projected along the a 
direction. Two different cation layers labeled A and B, respectively, are evidenced. Atomic ar- 
rangement of the ET molecules plane and band structure calculations relevant to (b) insulating 
layer A and (c) metallic layer B. (d) Fermi surface relevant to layer B. Dashed line marks the top 
of the highest occupied band and the Fermi level in (b) and (c), respectively. Labels Ao and Ai 
stand for the band extrema. F, X, Y, M and S refer to (0, 0), (a72, 0), (0, 6*72), {a* /2, b* /2) and 
(a*/2, -b*/2), respectively. 



where <l>o = h/e is the quantum flux and uq = 27T'^mekB /eh = 14.694 T/K. Effective 
masses are expressed in rrie units and chemical potential fi in Tesla. The index r] stands 
for all the possible classical orbits a, /3, /3 + a, 2/3 - a, etc., excluding their harmonics, 
with frequencies = m^(/i - A^) and effective masses m^. Frequencies depend on the 
chemical potential /i since they represent a filling factor corresponding to classical orbits. 
Coefficient Crj is the symmetry factor of the orbit rj. For example, Ca = C(3 = C2(3-a = 1 
and Cf3xa = 2. (pr] is the phase of the orbit r^, i.e. 7r/2 times the number of extrema of the 
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orbit in the kx or ky direction, for example, = = ^2(3-a = tt and (p(3-\-a = 27r. These 
phases are important for determining the sign of amphtudes. In particular, we define two 
fundamental orbits with frequencies, = mi(/i - Ai) {i.e. rua = ^i) and Ff^ = mo(/i - 
Ao) + mi(/i - Al) which can be written as Ff^ = m^{ii - A^) (see Fig. [TJd)). Therefore, 
mo + mi is identified to the effective mass m/3 of the MB orbit /3. The damping factors can 
be expressed as Rr^^p{T^B) = R^^pR^^pR^^p R^r^.p where the thermal. Dingle, MB and spin 
damping factors are given by: 



Bcos^sinh(-g^) 



R^'f = {ipoTHqor^ (4) 



respectively. Td is the Dingle temperature (Td = h/27rkBr where r is the relaxation 
time). Integers and are the (even) number of MB tunnelings and reflections, respec- 
tively, encountered by the quasiparticle along its closed trajectory. The MB tunneling and 
reflection probabilities are given by po = ex.p{—Bo/2B) and = 1 - Pq, respectively, where 
^0 is the MB fleld. g* is the effective Lande factoi0. The electron density per surface area 
He is given by dl^/d/i = -ng. In zero-fleld, Eq.[T] yields ng = (mo + mi)fio - ^oAq - miAi. 
In presence of a magnetic fleld, /i satisfles the following implicit equation: 



B cos 1 . F„ 

The oscillating part of the magnetization can be computed systematically by inserting 
Eq.[6]in Eq.[T]and deriving the resulting expression with respect to the magnetic fleld. Direct 
numerical resolution of this implicit equation could be done. Besides, a more user-friendly 
controlled expansion in term of damping factors Rr^^p can be derived systematically at any 
possible order. Keeping successive terms in powers of Rr^^p up to the second order, yields 
analytically the following expression for the oscillating part of the magnetization: 



5 cos (9 



77,77' p,p'^l 



sm I 27r — ^ + pcpr^ + p (prj> 

B cos 



sin 27r 



BcosO 



(7) 



where frequencies are evaluated at /i = /io, namely, = mrj{fio — ^r])- Eq. [7]can be 
written as a sum of periodic functions with frequencies Fi and Fourier amplitudes Ai where 
the index i stands for either classical orbits r] or 'forbidden' orbits such as /3 - a. Relevant 
examples for experimental analysis are given below: 

■•^FS warping can be accounted for by the additional damping factor Jo{4:7Tpm'qT±/ehB), where t± is the 
interlayer transfer integral [12 . In most q-2D organic metals, this factor is close to 1. 
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Fa Fa Fa 2 Fa 
Aa = Ra,l + Ra,lRa,2 + Ra,2Ra,3 H ^/3,l^a+/3,l , (8) 

^2q; = —- — - — Ra,2 H ~Ra 1 ~ n ~Ra,lRa,3^ (9) 

= -^i^/3,l + :7^i^/3,li^/3,2 + ^^i^/5,2i^/3,3 + ^i^/3,li^c.+/3,l. (10) 

^2/? = --^i?/5,2 + -^i?|,i-^^%l%3, (11) 

27rm/3 7rm/3 ^' 37rm/3 

= ^—^Ra,lRf3,l ^^a,2^a+/3,l ^^/3,2^q;+/3,1 , (12) 

7rm/3 7rm/3 7rm/3 

^^+q; = ^-^^/3+q;,1 H ^^^a, 1^/3,1, (13) 

^2/3-a — ^2/3-a,l + Ti ^a, 1^/3, 2 + 7 ^a,3^Q;+/3,2- U^j 

In the framework of the LK model, the amplitude of a Fourier component linked to a basic 
orbit involves only one damping factor [e.g. Aa oc Ra)- This is also the case for a MB orbit, 
the effective mass of which is given by the Falicov-Stachowiak model {e.g. m2(3-a = - 
^a) According to this statement, the leading term of Eqs. [8l and [TOl which accounts for 
the a and j3 components, respectively, reproduces the LK equation. Furthermore, deviations 
from the LK behaviour due to the high order terms involving products of damping factors 
are only significant in the low temperature and high field ranges. This statement also stands 
for 2/3-a since R2^-a is significantly higher than the product Ra,iR(3,2 (see Eq. [M]). As a 
consequence, the Falicov-Stachowiak model should apply for this component as well. In 
contrast, algebraic sum of damping terms in Eqs. [9land[T3l where minus signs accounts for 
TT dephasings, may cancel at field and temperature values dependent on the effective masses. 
Dingle temperature, MB field, etc. as displayed in Fig. [2] relevant to 2a. This point has 
already been reported for the second harmonic of the basic orbits both for compensated 
[13] and un-compensated [14] metals. In that respect, since algebraic sums are involved, 
care must be taken of the spin damping factor sigro. Indeed, the value of the spin damping 
factor, which is changed through the value of g* in Fig. [2l can have a drastic effect on 
the field and temperature dependence of the amplitude in this case. Finally, the Fourier 
amplitude of the 'forbidden orbit' /3 - a is accounted for by Eq. [121 Although a non-LK 
behaviour is evidenced, a monotonic field and temperature dependence is expected owing 
to the negligibly small value of the high order terms. The consistency of Eqs. [HI to [T^ with 
experimental data is examined in the next section. 

Results and discussion. — Even though the organic metal 6>-(ET)4CoBr4(C6H4Cl2) 
undergoes a phase transition at 343 K [6 1 , no additional phase transition is revealed by X-ray 
diffraction below room temperature, at which crystals are synthesized, down to 100 K. As 
displayed in Fig.[Tl the central feature of the crystalline structure at room temperature is that 
the unit cell has two different donors planes labeled A and B. According to band structure 
calculations, layer A is insulating while layer B is metallic. FS of layer B is composed of one 
closed tube centered at Y with an area of 18 % of the first Brillouin zone (FBZ) area (yielding 
a orbit in magnetic field) and two q-lD sheets with rather small gap in between. Therefore, 
it is expected that the MB orbit /3, which corresponds to the hole tube centered at F with 
an area equal to that of the FBZ is observed in moderate fields, giving rise to the textbook 



^ Since the LK equation involves a single damping term, the absolute value of the spin damping factor is 
considered in the case where data are analyzed on this basis. However, tt dephasing is observed on either 
side of spin-zero angles as reported in [TS) . 
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Fig. 2: (Colour on-line) Temperature dependence (mass plot) of the amplitude of the component 
2a. Symbols are data deduced from torque measurements of Fig. SJa) at 48.3 T. Solid lines are 
deduced from Eq.[9]with rria — 1.81, — 3.52, Td = 0.79 K and Bq — 35 T. Spin damping factor 
(see Eq. (5]) is either neglected {R^ = 1, p = 1 to 3) or accounted for by various g* values. All the 
data are normalized to their values at 12 K. 

linear chain of coupled orbits. In addition, since layers A are insulating, conducting planes 
are strongly separated from each other which suggests a strongly 2D behaviour. In line 
with this statement, the room temperature resistivity anisotropy p(interlayer)/p (in plane) 
is as large as 10^ [6]. Furthermore, according to the data in Fig. [3l this anisotropy further 
increases up to p(interlayer)/p(in plane) = 10^ as the temperature decreases down to liquid 
helium temperatures. 

Magnetic torque (dHvA oscillations) and tunnel diode oscillator [Shubnikov-de Haas 
(SdH) oscillations] data are displayed in Fig. HJ Besides Fourier components linked to the 
closed orbit a and to the MB orbit /3, few frequency combinations typical of the FS displayed 
in Fig. [T]are observed in the high field range of the torque spectra. In contrast, numerous 
combinations are observed in SdH spectra. As an example, the Fourier component 2/3 - 3a 
is clearly detected in Fig. ID^b). This frequency, which to our knowledge is observed for the 
first time, could correspond to a quantum interferometer path. However, its amplitude is 
unexpectedly large. The frequencies linked to the j3 and a orbits are Ff^iO = 0) = 4600 =b 
10 T, in agreement with the FBZ area at 100 K, and F^{e = 0) = 944 ± 4 T, respectively. 
The latter value yields a cross section area of 20.6 % of the FBZ area, in satisfactory 
agreement with the band structure calculations of Fig. [H Field- and temperature-dependent 
Fourier amplitudes can be studied at the light of the results given in the preceding section, 
keeping in mind that dHvA amplitudes A^ are related to torque oscillation amplitudes 
AL / BidiiiiO) . Effective masses, derived from Eq. [8] at various magnetic field 
values are field-independent: m{a) = 1.81 ± 0.05, m(/3) = 3.52 ± 0.19 (see Fig. [5fa)). In 
agreement with results of Ref. [16], m{/3) is twice m{a) within the error bars. At odds with 
torque results, effective masses deduced from SdH data, in particular those relevant to the 
P component, definitely increase as the magnetic field increases. This behaviour indicates 
that the above calculations, relevant to dHvA oscillations, cannot account for the SdH data 
at high field. 
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Fig. 3: (Colour on-line) Temperature dependence of zero-field in-plane (blue lower curve) and 
interlayer (black upper curve) resistance normalized to the values at 290 K. 

Focusing on dHvA oscillations, Dingle plots of a and f3 can be obtained plotting A^/ BR^^i^ 
where R^^ is given by Eq. O assuming the LK formalism accounts for the data, i.e. the 
high orders terms of Eq. [8] can be neglected. As reported in Fig. O all the data relevant 
to the a and (3 components lies, respectively, on the same curve within the experimental 
scattering which suggests that the above assumption holds. Dingle temperature and MB 
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F (kT) F (kT) 

Fig. 4: (Colour on-line) Fourier analysis in the field range 35-55.3 T of (a) the magnetic torque 
data (dHvA oscillations) and (b) the tunnel diode oscillator data (SdH oscillations) reported in the 
respective insets. Solid triangles are marks calculated with Fa(0 = 0) = 0.944 kT and F^{0 = 0) = 
4.60 kT. Data have been shifted down from each other by a constant amount. 



p-7 



A. Audouard et al. 




-8.5 



a 



b) 





'9.0 


' V 

1,4 k \ 


6 60 : 


} 


-9.5 


1,9 K 


k -S.65 - 






; iL.i i\ 


\-&.70 - 




10.0 


^ 3.5 K 
■ 4,2 K 






VB (T '} 



0.02 



0.03 0.04 
1/B (T') 



O.OS 



Fig. 5: (Colour on-line) (a) Effective masses vs. magnetic field for the a (circles) and (3 (squares) 
oscillations deduced from best fits of Eq.[8]to the torque (solid symbols) and tunnel diode oscillator 
(open symbols) data. Horizontal lines stand for the mean values of the effective masses deduced 
from torque data {rria, = 1.81, — 3.52). (b) Dingle plot of the a and /3 torque oscillations. 
where rf — (5^ is calculated with rria and mp deduced from (a). Solid lines are best fits of Eq. [8] 
and [10] to the data relevant to a and /3, respectively, obtained with Bo = 35 T and Td = 0.79 K 
(see Eqs. [2]to[4]). Inset displays an enlarged view of the fits for a. 



field can be derived for components which do not involve tunneling (n^ = in Eq. [5]). Al- 
though this is the case of a, a very broad maximum is observed in Fig. [5l leading to a very 
large uncertainty on the MB field. More specifically, Bq values ranging from to 70 T can 
account for the data, in which case Tjj ranges from 1.4 to 0.5 K, respectively. Assuming the 
Dingle temperature linked to the a orbit is the same as that of /3 considerably reduces the 
uncertainty since the well defined slope = -(iio^D^/3+25o) of the Dingle plot for /3 yield 
an additional relationship between T^^ and Bq. This assumption yields 5o = 35 ± 5 T and 
Td = 0.79 ± 0.10 K. With these values of effective masses. Dingle temperature and MB field, 
it is possible to check whether or not the high order terms of Eq. [8] can be neglected. Inset 
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Fig. 6: (Colour on-line) (a) Mass plots and (b) Dingle plots of /S - a (circles), /S -\- a (squares) and 
2/3 - a (diamonds). Solid lines are bests fits of Eqs. I12[ I13l and [Ml calculated with nia = 1.81, 
= 3.52, Bo = 35 T, Td = 0.79 K and ^* = 2. 
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of Fig.[5fb) displays the field dependence of A^/BR^^ calculated for the highest field range 
covered by the experiments, assuming further that g* = Although the contribution of 
these factors rises as the field increases, it remains smaller than the experimental scattering 
for both a and /3. 

The field and temperature dependence of the Fourier amplitude relevant to frequency 
combinations is displayed in Fig. [6l Since their amplitude is much smaller than that of the 
orbits a and /3, only mass plots at high field and Dingle analysis at the lowest explored 
temperature are considered. As already discussed for the second harmonics 2a (see Fig. [2]), 
amplitudes are strongly dependent on the spin damping factors value. Solid lines in Fig. [6] 
are the best fits of Eqs.[T2]to[T4j obtained with the same values of mo,, Tjj and as 

above, with only a prefactor as a free parameter for each of the data: a very good agreement 
is obtained, in particular for the 'forbidden orbit' /3 - a. It can be remarked that the 
behaviour of both 2/3 - a and P -\- a are close to the predictions of the Falicov-Stachowiak 
model. Whereas this feature is mainly due to the small value of Ra,iR^,2 compared to 
R2(3-a,i in the former case, it can be attributed to the peculiar set of values of the spin 
damping factors (^^,1^^, 2/^^+0;,! ~ -0-22) in the latter case. 

Summary and conclusion. — Analytical formulae for dHvA oscillations relevant to a 
class of FS illustrating the Pippard's model and observed in many q-2D organic conductors 
have been derived systematically taking into account the chemical potential oscillations in 
magnetic field (Eqs. [8] to [T4|) . Although high orders terms are observed, basic (a) and 
MB-induced {/3 and 2/3 - a) orbits can be accounted for by the LK and Falicov-Stachowiak 
semiclassical models, at least at moderate fields and temperatures. In contrast, the 'for- 
bidden orbit' P - a is accounted for by a non-LK equation, involving products of classical 
amplitudes relevant to a and p. Furthermore, strongly non-monotonic field and temperature 
dependence may be observed for the second harmonics of basic frequencies such as 2a and 
the MB orbit /3 + a, depending on the value of the spin damping factor Rf^. 

These formulae are in agreement with the dHvA oscillation spectra of the strongly 2D 
organic metal 6>-(ET)4CoBr4(C6ll4Cl2), assuming g* = 2. However, SdH oscillations spectra 
exhibit a tremendous number of Fourier components suggesting, besides quantum interfer- 
ence contributions, more severe deviations from the LK model that remains to be explained 
within a SdH theory including quantum interference contribution. 
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